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The  discontinuous  Galerkin  method  (DGM)  has  become,  in  recent  times,  one  of  the  most  widely  researched  and  utilized  discretization 
methodologies  for  solving  problems  in  science  and  engineering.  Fundamentally  based  upon  the  mathematical 

framework  of  variational  methods,  the  DG  methodology  provides  hope  that  computationally  fast,  efficient  and  robust  methods  can  be 
constructed  for  solving  real-world  problems.  Through  a  combination  of  a  dual  path  to  convergence  -  allowing  naturally  both  conforming  and 
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specific  solvers.  By  not  requiring  that  the  solution  be  continuous  across  element  boundaries,  the  DGM  provides  a  flexibility  that  can  be 
exploited  both  for  geometric  and  solution  adaptivity  and  for  parallelization. 

The  goal  of  this  proposal  is  to  research,  develop  and  implement  visualization  techniques  that  are  specifically  designed  to  respect  the 
mathematical  nature  of  high-order  DGM  simulation  fields.  As  our  intended  use  for  this  research  is  both  simulation  debugging  and  scientific 
exploratory  visualization,  we  seek  methods  and  implementations  that  act  on  high-order  DG  data  directly,  that  are  interactive,  and  that  have 
quantifiable  visual  error.  This  effort  is  designed  to  build  upon  our  previous  research  into  the  visualization  of  high-order  (continuous 
Galerkin)  finite  element  fields,  but  presents  a  distinct  challenge  in  that  we  must  now  consider  the  additional  computational  (e.g.  precision) 
and  visualization  challenges  of  element-wise  discontinuities. 
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Scientific  Progress 


Approach 


The  thesis  of  this  work  that  in  turn  drives  our  approach  is  that  a  visualization  methodology  for  high-order  finite  element  data  that 
exploits  the  high-order  nature  of  the  data  in  its  native  form  provides  a  visual  representation  that  introduces  no  (or  quantifiable) 
approximation  error  due  to  the  visualization  technique.  There  are  many  possible  reasons  why  this  thesis  is  beneficial  to  the 
high-order  finite  element  community.  The  proposed  visualization  techniques  (1)  use  the  data  in  its  native  form,  hence  helping  to 
allay  the  computational  scientists’  concern  that  information  is  being  lost  when  current  visualization  techniques  are  applied;  (2) 
allow  the  computational  scientists  to  focus  their  efforts  on  elimination  of  other  sources  of  error  (modeling  errors,  numerical 
(simulation)  errors, etc.)  because  it  eliminates  visualization  approximation  error;  and  (3)  provide  “ground  truth”  images  to  which 
low-order  visualization  techniques  applied  to  high-order  finite  element  data  can  be  compared.  To  facilitate  our  efforts,  we 
continue  to  build  upon  our  current  interactive,  nVIDIA  (www.nvidia.com)  Graphical  Processing  Unit  (GPU)  based  high-order 
FEM  visualization  software  called  EIVis(“Element  Visualizer”).  A  concrete  outcome  of  this  work,  in  addition  to  peer-reviewed 
publications  targeted  to  the  applied  mathematics,  visualization  and  application-domain  communities,  is  the  implementation  and 
dissemination  of  our  results  within  this  publicly  available  software  package. 

Scientific  Opportunities  and  Barriers 

Many  current  scientific  visualization  techniques  applied  to  higher-order  solutions  are  inadequate  when  used  for  knowledge 
extraction  and  assistance  in  reducing  the  error  budget  because  they  transform  high-order  data  to  low-order  representations  for 
visualization  purposes  -  a  process  which  in  and  of  itself  adds  “visualization  error”  to  the  error  budget.  The  scientist  is  currently 
burdened  with  determining  whether  or  not  an  anomaly  found  in  an  image  generated  by  a  visualization  technique  is  from  the 
modeling  and  discretization  assumptions  made  as  part  of  the  simulation  or  as  part  of  the  visualization  technique.  In  most  cases 
that  burden  is  high.  In  this  proposed  work  we  will  examine  several  common  scientific  visualization  techniques  used  for  depicting 
scalar  field  values  from  the  results  of  high-order  finite  element  (continuous  Galerkin  and  discontinuous  Galerkin)  simulations. 

We  propose  to  create  “highorder  cognizant”  visualization  algorithms  and  strategies  for  high-order  continuous  and  discontinuous 
data  on  planar  and  curved  elements. 

Significance 

The  proposed  research  impacts  three  areas:  the  mathematical  sciences,  the  computer  sciences,  and  the  interdisciplinary  bridge 
lying  between  these  two  areas.  The  high-order  finite  element  community  will  benefit  from  this  effort  through  the  proposed 
development  of  algorithms,  which  accurately  and  efficiently  render  simulation  results.  The  visualization  community  will  benefit 
through  the  exposure  to  the  high-order  finite  element  community  and,  in  particular,  the  numerical  methods  prevalently  found 
there.  Other  current  projects  funded  by  ARO’s  program  in  computational  mathematics  that  use  high  order  finite  elements  will 
find  immediate  utility  of  the  algorithms  and  implementations  discussed  herein  (e.g.  projects  funded  at  Brown  University,  RPI, 
etc.). 

Accomplishments 

This  report  summarizes  three  years  of  ARO  funding  on  the  current  project.  In  the  first  year,  a  graduate  student  (Mr.  Yichen 
Zhou)  began  working  on  the  project.  His  focus  during  the  first  year  was  learning  the  ElVis  software  and  beginning  to  examine 
mathematical  mechanisms  to  accelerate  ray-tracing  of  high-order  discontinuous  Galerkin  fields,  in  particular  whether  high-order 
root  finding  combined  with  interval  arithmetic  can  be  used  within  our  ElVis  GPU  environment.  Mr.  Zhou  produced  a  project 
report  for  his  efforts  and  gave  a  group  presentation  of  his  work.  After  assessment  by  both  Pis  of  his  work  and  performance  on 
the  project,  Mr.  Zhou  was  removed  from  the  project. 

To  bridge  our  efforts  while  searching  for  a  new  student,  we  engaged  two  post  doctoral  researchers,  one  on  the  Utah  side  and 
one  at  MIT.  On  the  MIT  side,  Dr.  Marshall  Galbraith  began  working  on  “ElVis  hardening”.  The  hardening  effort  consisted  of 
streamlining  the  compilation  procedure,  fixing  issues  relayed  to  us  by  MIT  users  (MIT  Project  X  users  who  provided  feedback 
during  the  evaluation  states  of  the  work),  working  on  in-code  commenting  and  on  documentation.  He  has  also  helped  to 
strategize  on  how  to  make  ElVis  more  cross  platform  stable.  His  efforts  with  respect  to  hardening  of  the  software  continued  until 
the  end  of  the  project,  and  led  to  a  second  release  of  the  software  which  was  better  documented  and  more  stable. 

On  the  Utah  side,  Dr.  Shankar  Sastry  worked  on  more  theoretical  components  of  the  ElVis  work  -  namely  questions  arising 
from  meshes.  The  original  question  posted  to  Dr.  Sastry  was  a  visualization  question:  could  we  accelerate  the  visualization  of 
high-order  methods  by  making  assumptions  on  the  inverse  mapping  used  between  world-space  (engineering  space)  and  the 
reference  space  (where  basis  functions  are  defined).  In  the  course  of  investigating  this  topic,  Dr.  Sastry  was  able  to  prove  some 
properties  about  quadratic  element  interpolation  that  did  not  exist  previously  in  the  literature.  A  paper  on  this  work  was 
published  at  the  International  Meshing  Roundtable  (IMR),  October  2013.  The  title,  author  list  and  abstract  are  as  follows: 

“On  Interpolation  Errors  over  Quadratic  Nodal  Triangular  Finite  Elements”  by  Shankar  P.  Sastry  and  Robert  M.  Kirby. 


Abstract:  Interpolation  techniques  are  used  to  estimate  function  values  and  their  derivatives  at  those  points  for  which  a 
numerical  solution  of  any  equation  is  not  explicitly  evaluated.  In  particular,  the  shape  functions  are  used  to  interpolate  a  solution 
(within  an  element)  of  a  partial  differential  equation  obtained  by  the  finite  element  method.  Mesh  generation  and  quality 
improvement  are  often  driven  by  the  objective  of  minimizing  the  bounds  on  the  error  of  the  interpolated  solution.  For  linear 
elements,  the  error  bounds  at  a 

point  have  been  derived  as  a  composite  function  of  the  shape  function  values  at  the  point  and  its  distance  from  the  element’s 
nodes.  We  extend  the  derivation  to  quadratic  triangular  elements  and  visualize  the  bounds  for  both  the  function  interpolant  and 
the  interpolant  of  its  derivative.  The  maximum  error  bound  for  a  function  interpolant  within  an  element  is  computed  using  the 
active  set 

method  for  constrained  optimization.  For  the  interpolant  of  the  derivative,  we  visually  observe  that  the  evaluation  of  the  bound  at 
the  corner  vertices  is  sufficient  to  find  the  maximum  bound  within  an  element.  We  characterize  the  bounds  and  develop  a  mesh 
quality  improvement  algorithm  that  optimizes  the  bounds  through  the  movement  (r-refinement)  of  both  the  corner  vertices  and 
edge  nodes  of  high-order  meshes. 

Dr.  Sastry’s  work  was  presented  at  IMR  and  sparked  new  work  that  was  presented  at  ICOSAHOM  2014.  During  the  third  year 
of  the  project,  we  extended  the  ICOSAHOM  work  into  a  journal  submission  (which  is  accepted  upon  revision)  entitled: 

“Interpolation  Error  Bounds  for  Curved,  High-Order  Finite  Elements  and  their  Implications  on  Adaptive  Mesh  Refinement”  by 
Shankar  P.  Sastry  and  Robert  M.  Kirby. 

Abstract:  Mesh  generation  and  adaptive  refinement  is  largely  driven  by  the  objective  of  minimizing  the  bounds  on  the 
interpolation  error  of  the  solution  of  the  partial  differential  equation  (PDE)  being  solved.  Thus,  the  characterization  and  analysis 
of  interpolation  error  bounds  for  curved,  high-order  finite  elements  is  often  desired  to  efficiently  obtain  the  solution  of  a  PDE 
through  the  finite  element  method  (FEM).  Although  the  order  of  convergence  of  the  interpolation  error  is  know  for  such 
elements,  a  quantitative  formulation  of  the  quality  of  an  element  is  still  not  known.  Moreover,  adaptive  mesh  refinement  (AMR) 
for  high-order  meshes  is  still  carried  out  through  Hessian-based  techniques  used  for  linear  elements.  We  derive  an  alternative 
formulation  for  the  interpoltion  error  bounds  for  curved,  high-order  element  that  uses  the  Taylor  series  expansion.  Our  error 
bounds  are  a  composition  function  of  the  shape  function  and  the  reference-  to  physical-space  mapping  functions.  This 
formulation  can  be  used  to  measure  the  quality  of  a  high-order  element  and  also  for  AMR.  We  find  that  the  AMR  for  pth-order 
FEM  should  be  driven  by  the  (p+1)th-order  partial  derivatives  of  the  function  for  which  the  PDE  is  being  solved.  Our  analysis 
can  also  be  extended  to  solve  for  anisotropic  functions.  Additionally,  we  note  that  our  error  bounds,  rather  than  the  Legesgue 
constant-based  bounds,  should  be  used  for  node  placement  in  a  high-order  element  if  certain  properties  of  the  function  being 
interpolated  are  known.  Our  numerical  experiments  indicate  that  the  error  bounds  derived  using  our  technique  is  asymptotically 
similar  to  the  actual  error,  l,e„  if  our  interpolation  error  bound  for  an  element  is  larger  than  it  is  for  other  elements,  the  actual 
error  is  also  larger  than  it  is  for  the  other  elements. 

In  addition,  on  the  Utah  side,  Mr.  Gene  Payne  (software  developer  here  at  the  SCI  Institute)  and  Ms.  Harshitha  Venkata 
(graduate  student)  worked  with  Dr.  Galbraith  on  the  hardening  and  deployment  of  the  code.  Mr.  Payne  led  the  architecting  work 
while  Ms.  Venkata  tested  the  software  on  various  use  cases  provided  along  with  the  software. 

Finally,  we  made  efforts  to  finalize  the  publication  of  three  papers  funded  under  previous  ARO  funding: 

Blake  Nelson,  Robert  M.  Kirby,  Robert  Haimes,  “GPU-Based  Volume  Visualization  From  High-Order  Finite  Element  Fields”, 
IEEE  Transactions  on  Visualization  and  Computer  Graphics,  Vol.  20,  No.  1,  pages  70-83,  2014. 

Tiago  Etiene,  Daniel  Jonsson,  Timo  Ropinski,  Carlos  Scheidegger,  Joao  Comba,  L.  Gustavo  Nonato,  Robert  M.  Kirby,  Anders 
Ynnerman  and  Claudio  T.  Silva,  “Verifying  Volume  Rendering  Using  Discretization  Error  Analysis”,  IEEE  Transactions  on 
Visualization  and  Computer  Graphics,  Vol.  20,  No.  1,  pages  140-154,  2014. 

Blake  Nelson,  Robert  M.  Kirby  and  Steven  Parker,  “Optimal  Expression  Evaluation  Through  the  Use  of  Expression  Templates 
When  Evaluating  Dense  Linear  Algebra  Operations”,  ACM  Transactions  on  Mathematical  Software,  Vol.  40,  Issue  3,  pages  21: 
1-21:21,2014. 

Collaborations  and  Leveraged  Funding 

•  Collaborated  with  Professor  Claudio  Silva,  Formally  of  School  of  Computing,  University  of  Utah,  and  now  at  the  Department  of 
Computer  Science,  NYU-Poly,  on  the  “Verifiable  Visualization”  work.  Mr.  Tiago  Etiene  was  the  graduate  student  (funded  under 
grants  directed  by  Prof.  Silva). 

•  Leveraged  Former  ARO  Funding  W91 1 NF-08-1-0517  (Program  Manager:  Dr.  Mike  Coyle)  efforts  that  resulted  in  the 
successful  completion  of  Dr.  Blake  Nelson  (PhD,  Computing,  University  of  Utah;  completed  Summer  2012).  Dr.  Nelson 
continued  to  finalize  research  papers  submitted  under  the  previous  grant. 

•  Leveraged  NSF  Funding  (IIS-0914564)  to  fund  a  software  developer  (Mr.  Gene  Payne)to  aid  in  the  “hardening”  of  ElVis  for 


both  public  distribution  as  well  as  future  graduate  student  usage  and  augmentation. 

Conclusions 

Visualization  is  often  employed  as  part  of  the  simulation  science  pipeline.  It  is  the  lens  through  which  scientists  often  examine 
their  data  for  deriving  new  science,  and  is  the  lens  used  to  view  modeling  and  discretization  interactions  within  their  simulations. 
As  such,  visualization  techniques  need  to  be  designed  not  only  to  elucidate  the  features  or  phenomena  of  interest  within  the 
data,  but  also  to  be  compatible  and  complementary  with  the  type  and  means  of  generating  the  data.  One  such  category  of 
simulation  data,  high-order  finite  element  methods  (also  known  as  spectral/hp  element  methods)  using  either  the  continuous 
Galerkin  or  discontinuous  Galerkin  formulation,  has  reached  a  level  of  sophistication  such  that  they  are  now  commonly  applied 
to  a  diverse  set  of  real-life  engineering  problems.  Visualizations  of  high-order  finite  element  results  that  do  not  respect  the  a 
priori  knowledge  of  how  the  data  were  produced  and  which  do  not  provide  a  quantification  of  the  visual  error  produced  may 
undermine  the  scientific  process  as  isolating  where  errors  and  assumptions  are  introduced  into  the  process  is  critical. 

Technology  Transfer 
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Objective 

The  discontinuous  Galerkin  method  (DGM)  has  become,  in  recent  times,  one  of  the  most  widely  researched  and  utilized 
discretization  methodologies  for  solving  problems  in  science  and  engineering.  Fundamentally  based  upon  the  mathematical 
framework  of  variational  methods,  the  DG  methodology  provides  hope  that  computationally  fast,  efficient  and  robust  methods  can  be 
constructed  for  solving  real-world  problems.  Through  a  combination  of  a  dual  path  to  convergence  -  allowing  naturally  both 
conforming  and  non-confonning  (hanging  node)  non-overlapping  discretizations  of  the  solution  domain  combined  with  (possibly 
non-uniform)  polynomial  enrichment  (also  known  as  p-refinement)  -  the  DG  methodology  provides  a  rich  mathematical  starting 
point  for  the  development  of  domain-specific  solvers.  By  not  requiring  that  the  solution  be  continuous  across  element  boundaries,  the 
DGM  provides  a  flexibility  that  can  be  exploited  both  for  geometric  and  solution  adaptivity  and  for  parallelization. 

The  goal  of  this  proposal  is  to  research,  develop  and  implement  visualization  techniques  that  are  specifically  designed  to  respect  the 
mathematical  nature  of  high-order  DGM  simulation  fields.  As  our  intended  use  for  this  research  is  both  simulation  debugging  and 
scientific  exploratory  visualization,  we  seek  methods  and  implementations  that  act  on  high-order  DG  data  directly,  that  are 
interactive,  and  that  have  quantifiable  visual  error.  This  effort  is  designed  to  build  upon  our  previous  research  into  the  visualization 
of  high-order  (continuous  Galerkin)  finite  element  fields,  but  presents  a  distinct  challenge  in  that  we  must  now  consider  the  additional 
computational  (e.g.  precision)  and  visualization  challenges  of  element-wise  discontinuities. 


Approach 

The  thesis  of  this  work  that  in  turn  drives  our  approach  is  that  a  visualization  methodology  for  high-order  finite  element  data  that 
exploits  the  high-order  nature  of  the  data  in  its  native  form  provides  a  visual  representation  that  introduces  no  (or  quantifiable) 
approximation  error  due  to  the  visualization  technique.  There  are  many  possible  reasons  why  this  thesis  is  beneficial  to  the  high-order 
finite  element  community.  The  proposed  visualization  techniques  (1)  use  the  data  in  its  native  form,  hence  helping  to  allay  the 
computational  scientists’  concern  that  information  is  being  lost  when  current  visualization  techniques  are  applied;  (2)  allow  the 
computational  scientists  to  focus  their  efforts  on  elimination  of  other  sources  of  error  (modeling  errors,  numerical  (simulation)  errors, 


etc.)  because  it  eliminates  visualization  approximation  error;  and  (3)  provide  “ground  truth”  images  to  which  low-order  visualization 
techniques  applied  to  high-order  finite  element  data  can  be  compared.  To  facilitate  our  efforts,  we  continue  to  build  upon  our  current 
interactive,  nVIDIA  (www.nvidia.com)  Graphical  Processing  Unit  (GPU)  based  high-order  FEM  visualization  software  called  ElVis 
(“Element  Visualizer”).  A  concrete  outcome  of  this  work,  in  addition  to  peer-reviewed  publications  targeted  to  the  applied 
mathematics,  visualization  and  application-domain  communities,  is  the  implementation  and  dissemination  of  our  results  within  this 
publicly  available  software  package. 


Scientific  Opportunities  and  Barriers 

Many  current  scientific  visualization  techniques  applied  to  higher-order  solutions  are  inadequate  when  used  for  knowledge  extraction 
and  assistance  in  reducing  the  error  budget  because  they  transform  high-order  data  to  low-order  representations  for  visualization 
purposes  -  a  process  which  in  and  of  itself  adds  “visualization  error”  to  the  error  budget.  The  scientist  is  currently  burdened  with 
determining  whether  or  not  an  anomaly  found  in  an  image  generated  by  a  visualization  technique  is  from  the  modeling  and 
discretization  assumptions  made  as  part  of  the  simulation  or  as  part  of  the  visualization  technique.  In  most  cases  that  burden  is  high. 
In  this  proposed  work  we  will  examine  several  common  scientific  visualization  techniques  used  for  depicting  scalar  field  values  from 
the  results  of  high-order  finite  element  (continuous  Galerkin  and  discontinuous  Galerkin)  simulations.  We  propose  to  create  “high- 
order  cognizant”  visualization  algorithms  and  strategies  for  high-order  continuous  and  discontinuous  data  on  planar  and  curved 
elements. 


Significance 

The  proposed  research  impacts  three  areas:  the  mathematical  sciences,  the  computer  sciences,  and  the  interdisciplinary  bridge  lying 
between  these  two  areas.  The  high-order  finite  element  community  will  benefit  from  this  effort  through  the  proposed  development  of 
algorithms,  which  accurately  and  efficiently  render  simulation  results.  The  visualization  community  will  benefit  through  the  exposure 
to  the  high-order  finite  element  community  and,  in  particular,  the  numerical  methods  prevalently  found  there.  Other  current  projects 
funded  by  ARO’s  program  in  computational  mathematics  that  use  high-order  finite  elements  will  find  immediate  utility  of  the 
algorithms  and  implementations  discussed  herein  (e.g.  projects  funded  at  Brown  University,  RPI,  etc.). 


Accomplishments 


This  report  summarizes  three  years  of  ARO  funding  on  the  current  project.  In  the  first  year,  a  graduate  student  (Mr.  Yichen  Zhou) 
began  working  on  the  project.  His  focus  during  the  first  year  was  learning  the  ElVis  software  and  beginning  to  examine 
mathematical  mechanisms  to  accelerate  ray-tracing  of  high-order  discontinuous  Galerkin  fields,  in  particular  whether  high-order  root 
finding  combined  with  interval  arithmetic  can  be  used  within  our  ElVis  GPU  environment.  Mr.  Zhou  produced  a  project  report  for 
his  efforts  and  gave  a  group  presentation  of  his  work.  After  assessment  by  both  Pis  of  his  work  and  performance  on  the  project,  Mr. 
Zhou  was  removed  from  the  project. 

To  bridge  our  efforts  while  searching  for  a  new  student,  we  engaged  two  post-doctoral  researchers,  one  on  the  Utah  side  and  one  at 
MIT.  On  the  MIT  side,  Dr.  Marshall  Galbraith  began  working  on  “ElVis  hardening”.  The  hardening  effort  consisted  of  streamlining 
the  compilation  procedure,  fixing  issues  relayed  to  us  by  MIT  users  (MIT  Project  X  users  who  provided  feedback  during  the 
evaluation  states  of  the  work),  working  on  in-code  commenting  and  on  documentation.  He  has  also  helped  to  strategize  on  how  to 
make  ElVis  more  cross-platform  stable.  His  efforts  with  respect  to  hardening  of  the  software  continued  until  the  end  of  the  project, 
and  led  to  a  second  release  of  the  software  which  was  better  documented  and  more  stable. 

On  the  Utah  side,  Dr.  Shankar  Sastry  worked  on  more  theoretical  components  of  the  ElVis  work  -  namely  questions  arising  from 
meshes.  The  original  question  posted  to  Dr.  Sastry  was  a  visualization  question:  could  we  accelerate  the  visualization  of  high-order 
methods  by  making  assumptions  on  the  inverse  mapping  used  between  world-space  (engineering  space)  and  the  reference  space 
(where  basis  functions  are  defined).  In  the  course  of  investigating  this  topic,  Dr.  Sastry  was  able  to  prove  some  properties  about 
quadratic  element  interpolation  that  did  not  exist  previously  in  the  literature.  A  paper  on  this  work  was  published  at  the  International 
Meshing  Roundtable  (IMR),  October  2013.  The  title,  author  list  and  abstract  are  as  follows: 

“On  Interpolation  Errors  over  Quadratic  Nodal  Triangular  Finite  Elements”  by  Shankar  P.  Sastry  and  Robert  M.  Kirby. 

Abstract:  Interpolation  techniques  are  used  to  estimate  function  values  and  their  derivatives  at  those  points  for  which  a  numerical 
solution  of  any  equation  is  not  explicitly  evaluated.  In  particular,  the  shape  functions  are  used  to  interpolate  a  solution  (within  an 
element)  of  a  partial  differential  equation  obtained  by  the  finite  element  method.  Mesh  generation  and  quality  improvement  are  often 
driven  by  the  objective  of  minimizing  the  bounds  on  the  error  of  the  interpolated  solution.  For  linear  elements,  the  error  bounds  at  a 
point  have  been  derived  as  a  composite  function  of  the  shape  function  values  at  the  point  and  its  distance  from  the  element’s  nodes. 
We  extend  the  derivation  to  quadratic  triangular  elements  and  visualize  the  bounds  for  both  the  function  interpolant  and  the 


interpolant  of  its  derivative.  The  maximum  error  bound  for  a  function  interpolant  within  an  element  is  computed  using  the  active  set 
method  for  constrained  optimization.  For  the  interpolant  of  the  derivative,  we  visually  observe  that  the  evaluation  of  the  bound  at  the 
corner  vertices  is  sufficient  to  find  the  maximum  bound  within  an  element.  We  characterize  the  bounds  and  develop  a  mesh  quality 
improvement  algorithm  that  optimizes  the  bounds  through  the  movement  (r-refinement)  of  both  the  comer  vertices  and  edge  nodes  of 
high-order  meshes. 

Dr.  Sastry’s  work  was  presented  at  IMR  and  sparked  new  work  that  was  presented  at  ICOSAHOM  2014.  During  the  third  year  of  the 
project,  we  extended  the  ICOSAHOM  work  into  a  journal  submission  (which  is  accepted  upon  revision)  entitled: 

“Interpolation  Error  Bounds  for  Curved,  High-Order  Finite  Elements  and  their  Implications  on  Adaptive  Mesh  Refinement”  by 
Shankar  P.  Sastry  and  Robert  M.  Kirby. 

Abstract:  Mesh  generation  and  adaptive  refinement  is  largely  driven  by  the  objective  of  minimizing  the  bounds  on  the  interpolation 
error  of  the  solution  of  the  partial  differential  equation  (PDE)  being  solved.  Thus,  the  characterization  and  analysis  of  interpolation 
error  bounds  for  curved,  high-order  finite  elements  is  often  desired  to  efficiently  obtain  the  solution  of  a  PDE  through  the  finite 
element  method  (FEM).  Although  the  order  of  convergence  of  the  interpolation  error  is  know  for  such  elements,  a  quantitative 
fonnulation  of  the  quality  of  an  element  is  still  not  known.  Moreover,  adaptive  mesh  refinement  (AMR)  for  high-order  meshes  is  still 
carried  out  through  Hessian-based  techniques  used  for  linear  elements.  We  derive  an  alternative  formulation  for  the  interpoltion  error 
bounds  for  curved,  high-order  element  that  uses  the  Taylor  series  expansion.  Our  error  bounds  are  a  composition  function  of  the 
shape  function  and  the  reference-  to  physical-space  mapping  functions.  This  fonnulation  can  be  used  to  measure  the  quality  of  a 
high-order  element  and  also  for  AMR.  We  find  that  the  AMR  for  pth-order  FEM  should  be  driven  by  the  (p+l)th-order  partial 
derivatives  of  the  function  for  which  the  PDE  is  being  solved.  Our  analysis  can  also  be  extended  to  solve  for  anisotropic  functions. 
Additionally,  we  note  that  our  error  bounds,  rather  than  the  Legesgue  constant-based  bounds,  should  be  used  for  node  placement  in  a 
high-order  element  if  certain  properties  of  the  function  being  interpolated  are  known.  Our  numerical  experiments  indicate  that  the 
error  bounds  derived  using  our  technique  is  asymptotically  similar  to  the  actual  error,  I,e„  if  our  interpolation  error  bound  for  an 
element  is  larger  than  it  is  for  other  elements,  the  actual  error  is  also  larger  than  it  is  for  the  other  elements. 

In  addition,  on  the  Utah  side,  Mr.  Gene  Payne  (software  developer  here  at  the  SCI  Institute)  and  Ms.  Harshitha  Venkata  (graduate 
student)  worked  with  Dr.  Galbraith  on  the  hardening  and  deployment  of  the  code.  Mr.  Payne  led  the  architecting  work  while  Ms. 
Venkata  tested  the  software  on  various  use  cases  provided  along  with  the  software. 

Finally,  we  made  efforts  to  finalize  the  publication  of  three  papers  funded  under  previous  ARO  funding: 


Blake  Nelson,  Robert  M.  Kirby,  Robert  Haimes,  “GPU-Based  Volume  Visualization  From  High-Order  Finite  Element  Fields”,  IEEE 
Transactions  on  Visualization  and  Computer  Graphics,  Vol.  20,  No.  1,  pages  70-83,  2014. 

Tiago  Etiene,  Daniel  Jonsson,  Timo  Ropinski,  Carlos  Scheidegger,  Joao  Comba,  L.  Gustavo  Nonato,  Robert  M.  Kirby,  Anders 
Ynnerman  and  Claudio  T.  Silva,  “Verifying  Volume  Rendering  Using  Discretization  Error  Analysis”,  IEEE  Transactions  on 
Visualization  and  Computer  Graphics,  Vol.  20,  No.  1,  pages  140-154,  2014. 

Blake  Nelson,  Robert  M.  Kirby  and  Steven  Parker,  “Optimal  Expression  Evaluation  Through  the  Use  of  Expression  Templates  When 
Evaluating  Dense  Linear  Algebra  Operations”,  ACM  Transactions  on  Mathematical  Software,  Vol.  40,  Issue  3,  pages  21:1-21:21, 
2014. 


Collaborations  and  Leveraged  Funding 

•  Collaborated  with  Professor  Claudio  Silva,  Formally  of  School  of  Computing,  University  of  Utah,  and  now  at  the  Department 

of  Computer  Science,  NYU-Poly,  on  the  “Verifiable  Visualization”  work.  Mr.  Tiago  Etiene  was  the  graduate  student  (funded 
under  grants  directed  by  Prof.  Silva). 

•  Leveraged  Former  ARO  Funding  W91  INF-08-1-0517  (Program  Manager:  Dr.  Mike  Coyle)  efforts  that  resulted  in  the 

successful  completion  of  Dr.  Blake  Nelson  (PhD,  Computing,  University  of  Utah;  completed  Summer  2012).  Dr.  Nelson 
continued  to  finalize  research  papers  submitted  under  the  previous  grant. 

•  Leveraged  NSF  Funding  (IIS-0914564)  to  fund  a  software  developer  (Mr.  Gene  Payne)  to  aid  in  the  “hardening”  of  ElVis  for 

both  public  distribution  as  well  as  future  graduate  student  usage  and  augmentation. 

Personnel  Funded 

•  Robert  M.  Kirby  (Utah)  -  PI 

•  Robert  Haimes  (MIT)  -  co-PI 

•  Mr.  Yichen  Zhou  (Utah)  -  graduate  student  (1st  year  of  project) 

•  Ms.  Harshitha  Venkata  -  graduate  student  (3ld  year  of  project) 

•  Dr.  Marshall  Galbraith  (MIT)  -  post-doc 

•  Dr.  Shankar  Sastry  (Utah)  -  post-doc 

•  Mr.  Gene  Payne  -  software  developer 


Conclusions 


Visualization  is  often  employed  as  part  of  the  simulation  science  pipeline.  It  is  the  lens  through  which  scientists  often  examine  their 
data  for  deriving  new  science,  and  is  the  lens  used  to  view  modeling  and  discretization  interactions  within  their  simulations.  As  such, 
visualization  techniques  need  to  be  designed  not  only  to  elucidate  the  features  or  phenomena  of  interest  within  the  data,  but  also  to  be 
compatible  and  complementary  with  the  type  and  means  of  generating  the  data.  One  such  category  of  simulation  data,  high-order 
finite  element  methods  (also  known  as  spectral//^?  element  methods)  using  either  the  continuous  Galerkin  or  discontinuous  Galerkin 
formulation,  has  reached  a  level  of  sophistication  such  that  they  are  now  commonly  applied  to  a  diverse  set  of  real-life  engineering 
problems.  Visualizations  of  high-order  finite  element  results  that  do  not  respect  the  a  priori  knowledge  of  how  the  data  were 
produced  and  which  do  not  provide  a  quantification  of  the  visual  error  produced  may  undermine  the  scientific  process  as  isolating 
where  errors  and  assumptions  are  introduced  into  the  process  is  critical. 


Technology  Transfer 

None  to  report  at  this  time. 


